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We investigate the non-Markovianity of continuous variable Gaussian quantum channels through the evolu¬ 
tion of an operational metrological quantifier, namely the Gaussian interferometric power, which captures the 
minimal precision that can be achieved using bipartite Gaussian probes in a black-box phase estimation setup, 
where the phase shift generator is a priori unknown. We observe that the monotonicity of the Gaussian interfer¬ 
ometric power under the action of local Gaussian quantum channels on the ancillary arm of the bipartite probes 
is a natural indicator of Markovian dynamics; consequently, its breakdown for specific maps can be used to 
construct a witness and an effective quantifier of non-Markovianity. In our work, we consider two paradigmatic 
Gaussian models, the damping master equation and the quantum Brownian motion, and identify analytically 
and numerically the parameter regimes that give rise to non-Markovian dynamics. We then quantify the degree 
of non-Markovianity of the channels in terms of Gaussian interferometric power, showing in particular that even 
nonentangled probes can be useful to witness non-Markovianity. This establishes an interesting link between 
the dynamics of bipartite continuous variable open systems and their potential for optical interferometry. The re¬ 
sults are an important supplement to the recent research on characterization of non-Markovianity in continuous 
variable systems. 

PACS numbers: 03.65.Yz, 03.65.Ta, 42.50.Lc 


I. INTRODUCTION 

The study of the dynamics of open quantum systems has 
received a lot of attention in recent research [1-3], especially 
with regards to the nature of the interaction between the sys¬ 
tem and its surrounding environment. This is a crucial as¬ 
pect in quantum information theory, where quantum resources 
irreversibly decohere under nonunitary evolutions described 
through quantum channels, and are rendered less useful for 
quantum protocols. In ideal terms, the dynamics of open sys¬ 
tems can be defined via a weak system-environment coupling 
and a long system relaxation time-scale, that entails a non- 
retrievable transfer of information from the system to the en¬ 
vironment. Such dynamics are modelled using master equa¬ 
tions of the Lindblad form [4-6], and are called Markovian 
[7]. However, a more realistic description of open systems dy¬ 
namics incorporates a stronger coupling with the environment 
and operates at shorter relaxation times that are comparable to 
the environment correlation time-scales: they are, therefore, 
inherently non-Markovian, and cannot be described by Lind- 
bladian completely positive semigroup maps [8-10]. Non- 
Markovianity is an essential feature in the open dynamics of 
more complex quantum systems, such as photosynthetic pig¬ 
ment protein complexes [11, 12], quantum dots in photonic- 
crystal cavities [13], and other biological [14, 15] and many- 
body strongly correlated systems [9]. Further, non-Markovian 
dynamics allow for an active backflow of information from 
the environment to the system [16-19], which has significant 
implications in quantum information science as observed by 
its role in quantum key distribution [20], metrology [21, 22], 
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quantum Darwinism [23], preservation of quantum correla¬ 
tions [24, 25], thermodynamical work extraction [26], and in 
designing enhanced quantum protocols [27-29]. 

One of the primary difficulty in analyzing non-Markovian 
dynamics is the complicated mathematical description of the 
dynamical quantum maps, as compared to Markovian pro¬ 
cesses [10, 30, 31]. This currently limits the analytical so¬ 
lution of non-Markovian maps to only a few quantum models. 
However, these models are of immense importance as they of¬ 
fer a realistic characterization for a variety of open systems, 
with a huge potential in practical applications [20-25, 27- 
29, 32-34]. These theoretical investigations have been de¬ 
veloped, in parallel, with a substantial increase in engineered 
environments that allow experimental modelling and control 
of non-Markovian dynamics [35-42]. 

More generally, in recent years, there have been several 
studies focused on the characterization and possible quan¬ 
tification of non-Markovianity in open quantum system dy¬ 
namics [10, 43]. The main focus of several works on non- 
Markovianity has been to capture the deviation from the dy¬ 
namical semigroup properties of the completely positive trace 
preserving (CPTP) quantum maps [8, 9]. In other words, 
a defining feature of Markovian dynamics is the divisibility 
of the CPTP map, and any departure from this property is 
a valid indicator of non-Markovianity [9]. Hence, the non¬ 
monotonic behaviour under non-divisible CPTP maps of a 
suitable quantum information-theoretic quantity, such as dis- 
tinguishability [16], entanglement [9], or quantum mutual in¬ 
formation [44], can be used to witness and quantify non- 
Markovianity. Other definitions have exploited the dynamical 
behavior of Fisher information [45], accessible information 
[46], local uncertainty [47], and interferometric power [48] 
(for other witnesses and quantifiers, see [49-52]). However, 
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the various measures and witnesses do not reliably embody 
necessary and sufficient conditions for all non-Markovian pro¬ 
cesses, and in general may not be compatible with each other 
[10, 30, 31, 43]. It is thus desirable to have a more concrete 
theoretical formalism for non-Markovianity in quantum pro¬ 
cesses, and further studies are required to fill the gap in the 
contemporary understanding of the subject. 

Though the majority of studies have sought to address the 
question of non-Markovianity in the dynamics of discrete 
level open quantum systems, such as qubits and spin chains, 
an analogous characterization for continuous variable (CV) 
systems remains much less developed, despite the fact that 
CV systems, and in particular Gaussian states thereof, con¬ 
stitute fundamental and highly controllable resources for a 
plethora of quantum information and communication proto¬ 
cols [53-55]. However, in recent times, there have been a few 
attempts to extend the theoretical formalism to characterize 
non-Markovianity in discrete level systems to the CV scenario 
[40, 56-59]. Most notably, approaches to capture the non¬ 
monotonicity of distinguishability [56] and non-divisibility of 
dynamical maps [40, 59] have been used to define important 
criteria to witness and measure non-Markovianity of Gaussian 
quantum dynamical maps, i.e., maps preserving the Gaussian- 
ity of their inputs. Attempts have also been made to character¬ 
ize non-Markovianity in terms of the universality of Gaussian 
dynamical maps [58] and volume of the physical Gaussian 
state space [57]. 

In this paper, we seek to define a measure of non- 
Markovianity for Gaussian quantum channels using an op¬ 
erational figure of merit defined in the context of quan¬ 
tum metrology, namely the Gaussian interferometric power 
(GIF) [60-62]. The GIF quantifies the guaranteed precision 
achieved using bipartite Gaussian probes in a black-box inter¬ 
ferometry setting, where the generator of the phase shift to be 
estimated is a priori unknown. The non-monotonic evolution 
of the interferometric power has been very recently shown, by 
some of us, to be useful to construct a valid witness and quan¬ 
tifier of non-Markovianity in discrete level open systems, effi¬ 
ciently computable in the case of single-qubit dynamics [48]. 
In the present work, we extend the formalism to investigate the 
specific conditions that characterize the parameter regimes re¬ 
lated to Markovian and non-Markovian dynamics in Gaussian 
channels. We consider two paradigmatic CV Gaussian maps 
based on the damping master equation [63] and the quantum 
Brownian motion [40, 64]. We observe that our characteri¬ 
zation, based on GIF, allows one to define a computable wit¬ 
ness of non-Markovianity for the considered Gaussian dynam¬ 
ics, consistent with the necessary and sufficient conditions for 
non-Markovian dynamics provided by non-divisibility of the 
Gaussian maps. We find that non-Markovianity can be effi¬ 
ciently witnessed even with nonentangled probes. Further, we 
investigate the optimal Gaussian probes that are needed to ob¬ 
tain a quantitative measure of non-Markovianity and analyze 
the scaling of this quantity with the mean energy of the probes, 
thus discussing the role of entanglement for enhanced sensi¬ 
tivity in the detection of non-Markovianity. We also comment 
on the robustness of our results with varying channels param¬ 
eters, including the response to finite bath temperatures. 


The organization of the paper is as follows. In Sec. II, 
we present a brief account on Gaussian states and dynamical 
maps, followed by a description of the GIF measure. Next, in 
Sec. Ill, we present the characterization of non-Markovianity 
for Gaussian channels and we illustrate the framework of this 
paper in two examples: the damping master equation with a 
single decay parameter and the quantum Brownian motion. 
We conclude with a discussion on the results in Sec. IV. 


11. GAUSSIAN STATES, GAUSSIAN CHANNELS, AND 
GAUSSIAN INTEREEROMETRIC POWER 


In this section, we briefly introduce the basic mathemati¬ 
cal formalism to describe Gaussian quantum states and dy¬ 
namical maps that are relevant to CV quantum information 
and the present investigation. In our work, we study non- 
Markovianity in CV systems by focusing on Gaussian states 
and Gaussian channels and by analyzing the evolution of the 
GIF. 


A. Gaussian states 


A CV system of two modes A and B (with annihilation op¬ 
erators a and b respectively), can be defined by the quadrature 
vector O = where - (a/i -i- a|)/ V2 and 

pk - {o-k - a\)l yfli, where k — A, B (assuming natural units, 
b = 1). The quadratures obey the canonical commutation re¬ 
lations [Oj, Ok] - i0.jk, with the two-mode symplectic form 




0 i\®' 

-1 oj 


( 1 ) 


A Gaussian state pab [55, 65-67] is represented by a Gaus¬ 
sian characteristic function in phase space, and is completely 
characterized by its first and second statistical moments of 
the quadrature vector, given respectively by the displacement 
vector 5ab - i^j) and the covariance matrix (Tab = (o'jk), 
where 6j = tr[pABOj] and ajk = ir[pAB{(Oj - dj), (Ok - di)]+] 
(with j,k — 1,...,4 for two modes), and (■,■]+ is the anti¬ 
commutator. With no loss of generality, one can set the 
first moments as null and, for all informational purposes, any 
Gaussian state can be completely determined by its covariance 
matrix. A bona fide condition satisfied by all physical Gaus¬ 
sian states is the Robertson-Schrodinger uncertainty relation, 
given by 


(Tab + — 0 - ( 2 ) 

Its worth recalling that, by local symplectic operations (equiv¬ 
alent to local changes of basis on the state), every two-mode 
covariance matrix can be transformed to a standard form, with 
diagonal 2x2 subblocks that can be written as: 


(Tab 


7^ /3j’ 


(3) 
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where a - diagja, a), ji - Amg[b, b], y = diag{c, dj, such 
that fl, ^7 > 1, c > |£/| > 0. The total mean number of excita¬ 
tions (proportional to the total mean energy) of a two-mode 
Gaussian state can be defined as: E = riA + ng - 2n, where 
fiA - (tr[Q:] -2)/4 and ng - (tr[/3] -2)/4 are the mean number 
of excitations in modes A and B respectively, and h denotes the 
mean number of excitations per mode. Throughout this paper 
we will always impose the physical assumption that any initial 
state be constrained to a finite mean energy. 


B. Gaussian channels 

In this work we study the open dynamics of Gaussian states 
through Gaussian channels, i.e. quantum channels that map 
Gaussian states into Gaussian states. If the final evolved 
state pAg{t) is Gaussian, its properties can similarly be stud¬ 
ied through its covariance matrix. The dynamical evolution 
of a two-mode covariance matrix cr subjected to a Gaussian 
channel can be written as [59, 68, 69]: 

cr(0) ^ (T(f) = X{t) (T{0)X{tf + Y{t), (4) 

where X,Y are real 4x4 matrices (for n-mode systems, 2n x 
2n matrices). For completely positive maps, the matrices must 
satisfy [59, 68,69]: 

Y + in- ixnx^ > 0. (5) 

If, as in the case considered in this paper, only mode A is 
subjected to a Gaussian channel, the evolved two-mode state 
is mapped to 

= [( V^l (OIa) © lB]^CrAB(0)[( V^i(OIa) © Ifi] 

+ A2(f)lA © Og, (6) 

where (t, 4£(0) is in the standard form Eq. (3), is the 2x2 
identity matrix acting on mode A, Og is a null 2x2 matrix 
acting on mode B, and finally Aj (f) and A 2 (f) are terms related 
to the specific dynamical evolution. Explicitly, one can write 
the evolved covariance matrix as follows, 

' aAi(f) -I- A2(0 0 c V^i(0 0 

0 aAi(f) + A2(f) 0 dy/Ai(t) 

cVARO 0 b 0 ■ 

, 0 d VAi(f) 0 b , 

(7) 


C. Gaussian Interferometric Power 

The GIF quantifies the ability of a two-mode Gaussian 
probe to estimate a local phase shift in a worst case scenario 
[61, 62], according to an operational setting generally known 
as black-box interferometry [60]. In the CV Gaussian sce¬ 
nario, the protocol to define the GIF can be summarized as 
follows: a two-mode Gaussian state pAg is prepared by two 
parties Alice and Bob as a probe for an interferometer; mode 


B (which we call the ancilla) enters a black-box in which it un¬ 
dergoes a unitary transformation where the pa¬ 

rameter (f) is unknown, and only the spectrum of the generator 
Hg (assumed nondegenerate to avoid trivial dynamics) is ini¬ 
tially known. Fixing the spectrum of Hg to be harmonic, and 
restricting to Gaussianity preserving transformations, the uni¬ 
tary can be written as = Vg^gVE, where Wg - 
is a conventional phase shift operator and Vg is an arbitrary 
Gaussian unitary transformation. The transformed two-mode 
state is given by 


= (I^ ® U^g) PAg (Ia ® U^g^ . (8) 


The information on the black-box generator is provided to the 
two parties only after the transformation (i.e., the choice of Vg 
is disclosed), thus allowing for optimal measurements to be 
performed on the output. For any given setting of Vg, the ob¬ 
jective of the interferometric setup is to deduce the unknown 
phase (j) with the maximum possible precision, i.e., to con¬ 
struct the best possible estimator of the parameter. Assuming 
a large number k of copies of the probing state pAg are initially 
prepared, and letting the interferometric trial be repeated ac¬ 
cordingly K times, by collective processing one can construct 
an estimator (pest whose variance - {{pest - P)^) is con¬ 
strained by the Cramer-Rao bound. 


kAP^ > -;-, 

ripfgy 


(9) 


where 'FiPAg’^) is the quantum Fisher information, which can 
be defined as 


r(p>%) = -2 hm 

d<t>s0 


.^2 72/^0.t's Jt+dtpAss 

O ^V>AB ^PaB ) 

d{dpy 


( 10 ) 


with Fifing, p‘^g'^) being the Uhlmann fidelity. Since the 
inequality in (9) is asymptotically attainable, the quantum 
Fisher information directly quantifies the precision in the esti¬ 
mation of p for each given choice of the black-box generator, 
being inversely related to the minimum variance of the opti¬ 
mal estimator, and being by construction geometrically inter¬ 
preted as the rate (speed) at which the probing state changes 
following an infinitesimal change of the parameter to be esti¬ 
mated. 

In the black-box paradigm, the natural figure of merit one 
adopts to quantify the guaranteed precision allowed by a given 
probing state pab, in the estimation of a parameter p with in¬ 
complete prior information on the local generator is therefore 
given by a worst-case scenario, whereby one minimizes the 
quantum Fisher information over all the local generators (with 
the fixed spectrum). The resulting quantity is known as the in¬ 
terferometric power of pAB [60]. In the CV Gaussian setting, 
of relevance in this paper, the GIF of the two-mode Gaussian 
state Pab, with respect to mode B, is thus defined as: 


Q^gipAB)-Unf r(pt{h- (11) 

^ Vb 
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It has been proven that, both in finite-dimensional and in CV 
systems, the interferometric power is a measure of discord- 
type correlations for general mixed states reducing to a 
measure of entanglement in the particular case of pure states 
[60-62]. In the Gaussian case, the GIP satisfies in partic¬ 
ular the following properties [61, 62]: (i) it vanishes iff pab is 
a product state (since the only classically correlated Gaussian 
states are product states [70, 71 ]); (ii) it is invariant under local 
unitaries; (iii) it is monotonically nonincreasing under local 
CPTP operations on party A. The latter property is the cru¬ 
cial one which will allow us to construct a non-Markovianity 
indicator based on the GIP in the next Section. 

A handy feature of the GIP is that the minimization in 
Eq. (11) can be solved exactly for two-mode Gaussian states 
and local Gaussian unitaries. The result is a compact formula 
to calculate the GIP in terms of the covariance matrix ctab of 
an arbitrary two-mode Gaussian state, given by [61] 


Q^((^ab) 


Cx + + CyCi 

2C; 


( 12 ) 


where 


Cx - (h + h)(^ + Ii +-14) - il 
Cy - (I4 - 1)(1 -H /i +12 + 2/3 + I4) 

Q = ih + hXhh - h) + /3(1 + h){2h + h), (13) 

with the symplectic invariants of the general covariance ma¬ 
trix (3) given by I\ - deta, I 2 - det^. It, - det 7 , and 
I4 — dtlCTAB- 


III. CHARACTERIZATION OF NON-MARKOVIANITY IN 
GAUSSIAN CHANNELS 

A. Non-divisibility of Gaussian channels 

In this section, we focus on the non-Markovian character¬ 
istics of a Gaussian quantum channel in the open system dy¬ 
namics of CV states. The time evolution of an initial Gaussian 
state p(0), under a Gaussian dynamical map A,, is given by 
the relation p{i) - A,p{0), where A, is CPTP. Under Marko¬ 
vian dynamics. A, can be described by Lindbladian dynam¬ 
ical semigroups and the quantum map is divisible. Hence, 
a Markovian channel can be written as A, = A, ,/A,/ 0 , for 
t > t' > 0 [10, 31], which implies that the CPTP map describ¬ 
ing the dynamics of the system from time 0 to f can be divided 
into two CPTP maps from initial time 0 to some intermediate 
time f', and from f to the final time t. This implies that any 
quantum informational quantity that is by definition contrac¬ 
tive under CPTP maps, must be strictly monotonic under a 
Markovian evolution. 

Hence, either by observing a violation of the divisibility 
criterion, or by witnessing a breakdown of monotonicity of a 
contractive quantum information measure, one can get natural 
witnesses of non-Markovianity; a few proposals along these 
lines have been reported recently for CV systems [40, 56-59]. 
In our work, we use the non-monotonic evolution of the GIP 



Figure 1. (Color online) A schematic illustrating the introduced pro¬ 
tocol to characterize continuous variable non-Markovianity via the 
Gaussian interferometric power of a two-mode Gaussian state. The 
system mode A is subjected to a local quantum Gaussian dynami¬ 
cal map, modelling the interaction with the reservoir, whereas the 
ancillary mode B is subjected to an unknown Gaussian unitary trans¬ 
formation. The observation of a non-monotonic time evolution of the 
Gaussian interferometric power of the output two-mode state yields 
a valid witness of non-Markovian dynamics on the system mode. 


to detect and quantify Gaussian non-Markovianity. Before in¬ 
troducing our non-Markovianity indicator, we briefly recall a 
characterization of non-divisible Gaussian channels. 

Given a Gaussian dynamical map (see Sec. II), acting on the 
covariance matrix of a Gaussian state as described in Eq. (4), 
and defined by the matrices X and Y respecting the condi¬ 
tion of Eq. (5), it can be shown that a necessary and sufficient 
condition for non-Markovianity, due to a violation of the di¬ 
visibility criterion [9], is given by [59] 

Y(t + e) + ifl - iX(t + e)n X'^it + e)<0, (14) 


where e is any intermediate time instant and fl is the symplec¬ 
tic matrix defined earlier. While Eq. (14) is useful to detect a 
non-Markovian evolution, the authors of Ref. [59] have fur¬ 
ther proposed to quantify the degree of non-Markovianity in 
CV Gaussian channels by a measure defined as 


Nd = 




G(t)dt, 


(15) 


where the integral is performed over the time interval I, and 
G{t) = 2*2 limf^o+[|v*(f -H e, t)\ - Vk(t + e, f)], with v*(f -H e, t) 
being the eigenvalues of the symmetric matrix in the left hand 
side of Eq. (14). 


B. Witnessing and quantifying non-Markovianity via the 
Gaussian interferometric power 

We define an alternative measure of non-Markovianity in 
CV systems by making use of the non-monotonic behavior 
of the GIP under nondivisible CPTP Gaussian maps. The in¬ 
terferometric power has been shown to be a practical quan¬ 
tifier of non-Markovianity in the context of discrete variable 
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systems [48]. Here we follow a similar approach to investi¬ 
gate non-Markovianity in Gaussian CV systems; see Fig. 1 
for an illustration of the setup. Let us consider a Gaussian 
state pab of the two modes A and B, where the system mode 
A undergoes an open evolution (due to the presence of a reser¬ 
voir) modelled as a Gaussian CPTP map, while the ancillary 
mode B is sent through a black-box unitary. As described 
in Section IIC, the GIP for the output two-mode state Pab( 0 
can be exactly computed; this quantity captures the ability of 
such two-mode state to serve as a probe for the estimation 
of a phase shift imprinted during the black-box transforma¬ 
tion on the ancilla mode, in a worst-case scenario. Due to 
its contractivity under local channels on the system mode A, 
the GIP Q^ipABiO) is by definition a monotonically nonin¬ 
creasing function of time under any local divisible CPTP map 
acting on mode A. Hence, under Markovian Gaussian chan¬ 
nels, the time-derivative of the GIP is strictly non-positive, 
i.e., j^CPipABiO) ^ 0, for all f > 0. The non-Markovian char¬ 
acter in the dynamics of mode A is thus characterized by the 
violation of the monotonic behavior of the GIP. If one defines 
the quantity 


D{t) = j^QP(PAB(t)), (16) 

any positive value of D{t) captures the non-monotonic be¬ 
havior of the GIP under the considered Gaussian channel. 
Hence, for a given initial Gaussian state Pab(O) with a co- 
variance matrix a-ABiO), and denoting the evolved state as 
PABiO = (Aa ® Ib)p.4b( 0) following the action of a local Gaus¬ 
sian channel A^ on the system mode A, we define a metrolog¬ 
ical witness of Gaussian non-Markovianity by the quantity 

N^(A)^ f D{t)dt, (17) 

J®(0>0 

where the integral is over all the time instants t such that 
G(t) > 0. The superscript cr in Nq emphasizes the fact that 
the witness is specific to a given initial Gaussian state with a 
covariance matrix cr = cr^B(O). Following [9, 16], the non- 
Markovianity witness in Eq. (17) can be optimized over the 
set of all initial Gaussian states to obtain a proper measure of 
Gaussian non-Markovianity, defined as 

Nq(A)^ max7V^(A). (18) 

O'Afi(O) 

Numerically, the integral may be evaluated in small interval 
integrals: 


7Vo(A) = max 

O-AS(O) 


2 r 


Dk{t)dt. 




(19) 


In the following subsections we show how the above wit¬ 
ness and measure of CV non-Markovianity can be applied to 
investigate the characteristics of physical quantum maps. For 
the purposes of our investigation, we consider two fundamen¬ 
tal Gaussian channels that are based on the damping master 
equation [63] and the quantum Brownian motion [40, 64]. 
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Figure 2. (Color online) Measure of non-Markovianity Nq for the 
damping master equation, as a function of the mean number of ex¬ 
citations h and of the coupling constant a. All the quantities plotted 
are dimensionless. 



C. Non-Markovianity in the damping model with a single 
decay parameter 

We begin with a Gaussian channel characterized by a single 
time-dependent decay parameter, which is represented by a 
well-described damping master equation, given by 

“T - Q;-^^(2apa^ - (a^a,p]+), (20) 

dt 2 

where a is the coupling constant, and y{f) the time-dependent 
damping coefficient. It can be easily shown that for the damp¬ 
ing master equation the divisibility property is satisfied by the 
condition y{t) >0, V f > 0 [56, 59]. Thus, the damping mas¬ 
ter equation provides an excellent model to test the validity of 
any prospective non-Markovian quantifier. 

In our setting, the covariance matrix of a two-mode Gaus¬ 
sian state whose mode A is subject to the damping master 
equation of Eq. (20), is mapped as: 

ctabU) = [(e-"«/2i^)©lB]V^B(0)[(e-"«/2j^)©lB] 

-^(1-e-^<b)i^©OB, (21) 

where 

x{f) - a f 2y(s)ds. (22) 

Jo 

If the initial covariance matrix (t,4b( 0) is in standard form, as 
in Eq. (3), then after local damping the covariance matrix is 
mapped by Eq. (21) to the evolved state ccABit) given explicitly 

by 


CTABit) = 


a(t) 0 c(t) 0 

0 a(t) 0 d(t) 

c(t) 0 b 0 

0 d{t) 0 b , 


(23) 
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Figure 3. (Color online) Witness of non-Markovianity for 10^ ran¬ 
domly generated two-mode Gaussian states, cr,4B(0), captured by Nq 
as defined in Eq. (25), for the damping master equation. Nq is plot¬ 
ted as a function of the mean number of excitations it for a fixed 
a = 0.1. The (red online) dashed line represents the pure two-mode 
squeezed states and the (blue online) solid line is for the two-mode 
mixed thermal states. It is evident from the figure that pure two¬ 
mode squeezed states optimize the witness Nq, to realize the non- 
Markovianity measure Nq plotted in Fig. 2. All the quantities plotted 
are dimensionless. 


where a(t) = H- (1 - e'^®), c(t) = d{t) = 

q-Q illustrate quantitatively our characterization of 
non-Markovianity as measured using the GIF, we can choose 
on a specific instance of damping regime, given by [56] 



e sinf, 

e-'r/4^ 


if f < 5n/2 
if f > 57t/2. 


(24) 


In this case, there exists only a single interval, n < t < 2n, 
where y{t) is negative. Therefore, the non-Markovian witness 
Nq and measure Nq, based on GIF, can be analytically eval¬ 
uated by the expressions 

7V^(A) = Q%{aAB{t = 2n)) - Q%{cTAB{t = tt)), (25) 

Nq(N) = max N^. (26) 

<tas(0) ^ 

Figure 2 shows the measure of non-Markovianity, Nq, for 
the damping master equation as a function of the mean num¬ 
ber of excitations per mode, h, and of the coupling constant 
a. We note that, for fixed a, the value of Nq increases with 
higher values of the mean energy. Similarly, for a constant 
mean energy h, Nq increases monotonically with a. 

In the next subsection we provide more details on the op¬ 
timal input states to witness non-Markovianity, and we care¬ 
fully discuss the role of entanglement in the detection. 


D. Non-Markovianity witness beyond entangled probes 

A remarkable aspect of characterizing non-Markovianity 
using GIF is its ability to witness the non-Markovian dy¬ 
namics of a local Gaussian channel by using two-mode 
probes which exhibit quantum correlations beyond entangle¬ 
ment [61]. To illustrate this, we can focus on two important 


classes of two-mode Gaussian states: the mixed thermal states 
(MTS), which are always unentangled, and the squeezed ther¬ 
mal states (STS), which can be entangled or separable depend¬ 
ing on the trade-off between squeezing and noise. Both these 
classes of states can be easily engineered with the current tool¬ 
box of quantum optics, using e.g. single-mode squeezing and 
and linear optical operations, such as beam-splitters [55, 72]. 
In particular, the MTS can be generated by letting two single¬ 
mode thermal states, at different equilibrium temperatures, in¬ 
teract via a balanced 50:50 beam-splitter. The STS can be in¬ 
stead engineered by first squeezing two single-mode thermal 
states in complementary quadratures, and then letting them 
interact via a 50:50 beam-splitter. The covariance matrices of 
these two classes of states can be explicitly written as. 


^ A D 


k\e^''^x 0 

0 

0 

0 kie^'''y 


0 

0 

kie^''G 0 

0 k\e^''N 


and (27) 


k 2 x’ 

0 

k2y' 

0 

0 

k 2 x' 

0 

-hy' 

^ 2 / 

0 

k 2 x' 

0 

0 

-hy' 

0 

k 2 x' 


(28) 


for MTS and STS, respectively, where x = cosh(2ri), y - 
sinh(2ri), x' - cosh(2r2), andy' = sinh(2r2). The parameters 
in Eqs. (27) and (28) are related to the mixedness of the initial 
thermal states {ki,k 2 > 1) and the strength of the Gaussian op¬ 
erations (ri, r 2 > 0) such as squeezing (in the STS case) and 
beam-splitter interaction. Although the MTS states of Eq. (27) 
are always separable, and the STS of Eq. (28) are separable 
when k 2 {x' - y') > 1, both classes of states always contain 
more general forms of quantum correlations such as quantum 
discord and have in general a nonzero GIF [61, 62, 73, 74], as 
soon as they are not completely uncorrelated (i.e., as soon as 
y,y' + 0). Since the resource ensuring a precision in black¬ 
box interferometry is indeed associated with discord-type cor¬ 
relations, as quantified by the GIF, and not with entanglement, 
and since we are exploiting this figure of merit to investigate 
non-Markovianity, both classes of states (as well as any other 
non-product Gaussian state) are in principle able to provide 
a useful witness of non-Markovianity in generic local Gaus¬ 
sian channels. There is furthermore a computational advan¬ 
tage in focusing on the two considered classes of states, since 
for both of them the covariance matrix is in standard form, 
Eq. (3), with d - +c; in this case, the expression (12) for the 
GIF reduces to the much simplified form [61, 62] 


^b(<^Ab)I i/=rc — 


2(ah — N ± 1) 


(29) 


In the specific case of the two classes of states defined above, 
we have a - b - kie^’’' cosh(2ri) and c - kie^’’' sinh(2ri) for 
MTS, and a - b - k 2 cosh(2r2),c = ^2 sinh(2r2) for STS. As 
mentioned in Sec. II, to issue a fair comparison and to find 
an appropriate balance between the choice of parameters be¬ 
tween these two classes of Gaussian states, we will require the 
physical assumption that all states are parametrized in terms 
of finite initial mean energy E - 2h. 
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■ Squeezed vacuum states c Mixed thermal states 


Figure 4. (Color online) Witness of non-Markovianity given by Nq 
for the damping master equation, as a function of the mean number of 
excitations h and the coupling constant a, for the classes of separable 
two-mode mixed thermal states (blue online lowermost surface) and 
entangled two-mode squeezed states (red online uppermost surface), 
given respectively by Eqs. (27) and (28), with ki = k 2 = 1. All the 
quantities plotted are dimensionless. 


For the damping master equation, we first observe that, 
within the classes of STS and MTS chosen as initial probe 
states, the non-Markovianity witness Nq is analytically opti¬ 
mized (for any fixed h), when the states are as pure as pos¬ 
sible, i.e., for = ^2 = 1. In this limit, the STS reduce 
to pure entangled two-mode squeezed vacuum states (also 
known as twin-beams), with all the energy invested in squeez¬ 
ing, while the separable MTS maintain a nonzero degree of 
mixedness. We have performed an extensive numerical in¬ 
vestigation of the witness Nq over completely general ran¬ 
domly generated two-mode Gaussian input covariance matri¬ 
ces, in order to identify the globally optimal probes. From 
the numerical analysis presented in Fig. 3 (in the instance of 
a - 0.1), corroborated by the analytical optimization within 
the classes of STS and MTS states, we conclude that two¬ 
mode pure STS (i.e., squeezed vacuum states), spanning the 
(red online) dashed curve, are the globally optimal probes 
which attain the maximum in Eq. (18), thus achieving the 
quantification of non-Markovianity of the damping channel 
through the measure Nq. This has been verified numerically 
over an extended range of values of a and h, and the result ob¬ 
tained by using squeezed vacuum probes matches exactly with 
the evaluation of the measure in Fig. 2. Remarkably, though, 
we can see that the separable MTS states with ^2 = 1, span¬ 
ning the blue solid curve in Fig. 3, provide a suboptimal but 
still significantly nonzero witness of non-Markovianity for the 
damping master equation. It is interesting to observe that sev¬ 
eral randomly generated states, even entangled, were found to 
perform worse than the MTS probes, resulting in significantly 
lower values of the witness Nq. This shows quite importantly 
that a robust non-Markovianity witness for Gaussian channels 


can be generated (at least for the single decay damping model) 
using separable mixtures of thermal states, which are signif¬ 
icantly more economical to prepare than entangled states re¬ 
quiring high degrees of squeezing. 

One can then wonder, if entanglement is not needed to 
detect non-Markovianity of Gaussian channels in our setup, 
what potential benefit it may bring other than a quantitative 
enhancement of the witness. The answer comes directly from 
the metrological nature [60, 61] of our non-Markovianity in¬ 
dicator. In optical interferometry [75], a linear scaling of the 
quantum Fisher information with the mean energy per mode 
h denotes the shot noise, or standard quantum limit; separa¬ 
ble probes can never surpass this limit. Conversely, entangled 
probes can allow for an enhanced estimation, namely a scaling 
of the quantum Fisher information proportional to fp-, which 
is the so-called Heisenberg limit [75, 76]. Given that our wit¬ 
ness Nq is constructed directly from the quantum Fisher in¬ 
formation, we may adopt the same metrological terminology 
in order to characterize the efficiency of our non-Markovianity 
witness and measure. We may then observe that the entangled 
STS allow for non-Markovianity detection at the Heisenberg 
limit, leading to a non-Markovianity measure scaling quadrat- 
ically with the mean energy of the probes, while the separa¬ 
ble MTS (and any other separable state) can still detect non- 
Markovianity, but providing a witness scaling at most linearly 
with the mean energy. We have compared the witnesses de¬ 
rived from optimized STS and MTS in the damping master 
equation in Fig. 4 as a function of both a and the mean energy 
per mode h, generalizing the curves reported in Fig. 3. The 
quadratic versus linear scaling with h is evident in the com¬ 
parison, although both quantities are monotonically increas¬ 
ing functions of a and are nonzero as soon as a > 0. 


E. Non-Markovianity in the quantum Brownian motion 


In this section we consider the non-Markovian dynamics of 
Gaussian states under the master equation for quantum Brow¬ 
nian motion (QBM). Such a dynamics has been recently ob¬ 
served experimentally in an opto-mechanical setup [40]. Us¬ 
ing the interaction picture, with a secular and weak coupling 
approximation, the master equation for QBM acting on the 
system mode A is given by [64], 


dp 

dt 


a - (lapa' - [a'a,p]+) 


+Q' 


--- {2a^pa - {aa\p]+). 


(30) 


where a is the coupling constant, and the coefficients A(f) and 
7 (f) are the diffusion and damping term, respectively. The 
master equation for QBM is an approximation of the exact 
master equation [56, 77], 


dp 

dt 


= -A(f)[§A, [qA,p]\ + n(f)[§A, [pA,pW + 


+ y^it)[q\,p] - iy(t)[pA,[pA,p]+\, 


(31) 


where the time-dependent coefficients appearing in the equa¬ 
tion depend on the the reservoir spectral density, and are in- 
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Figure 5. (Color online) Non-Markovianity indicators in the quantum Brownian motion, plotted as a function of a, for different (scaled) 
reservoir temperatures: (a) T = 0, (b) T = 0.5, (c) T = 1, and (d) T = 4. The non-Markovianity witness JVg is plotted for the squeezed 
thermal state with fe = 1 (red online dashed lines) and the mixed thermal state with ki = \ (blue online solid lines). All the states have fixed 
average number of excitations, h = 2.5. The Ohmic cut-off frequency and the characteristic frequency have been set at ojc = 1 and ojg = 4, 
respectively. The measure based on divisibility criteria, Afo, is plotted as well for comparison (black online dotted lines). All the quantities 
plotted are dimensionless. 


terpreted as follows: R(t) phase shift; y{t) —> damping; 
A(f) (n(f)) —> normal (anomalous) diffusion. The diffusion 
and damping coefficients in the master equation given by 
Eq. (30) can be written as 

A(f)- I ds J" dajJ(a))^N(a)) +-'jcos(cL)os)cos(oJs); 

(32) 

y(t) - I ds I dtoJito) sin(wos) sin(a)s), (33) 
Jo Jo 

where N(a)) = {exp[fia)/kB'r] - 1)“^ is the mean number of 
thermal photons with frequency o), and cjq is the characteristic 
frequency of the system. We note that the coefficients A(f) and 
y{t) can be derived once the spectral density J{(jj) and temper¬ 
ature T of the reservoir are known. In our investigation, we 
consider the specific case of an Ohmic spectral density of the 
reservoir with an exponential cut-off frequency a>c, such that 
J(to) = 

Considering as usual an initial two-mode covariance matrix 
(Tab{0), and letting mode A undergo a QBM evolution under 
the master equation given by Eq. (30), the evolved covariance 
matrix is obtained as follows: 

f e-’^^^^A(s)ds(lAeOB), (34) 

Jo 


where x(t) - a 2y(s)ds, and the expressions of A(f) and y(t) 
are given by Eqs. (32) and (33), respectively. The dynamics 
of the Gaussian state under local QBM, given via Eq. (34), 
is mathematically more intricate, as compared to the case of 
Eq. (21) for the local damping channel, due to the presence of 
two decay components, A(f) and y(t). The non-Markovian dy¬ 
namics of the Gaussian system is determined by the interplay 
between these two terms. It is known that, for the QBM mas¬ 
ter equation, the non-divisibility condition given by Eq. (14) 
is satisfied for A(f) < \yit)\, or equivalently for negative values 
of and [56, 59]. 

Eor the specific Ohmic spectral density under considera¬ 
tion, the non-Markovianity in the QBM depends on the ra¬ 
tio (jOqIcjJc between the characteristic frequency of the system. 
Wo, and the cut-off frequency, Wc- In particular, one can ex¬ 
pect that in the regime wq/Wc ^ 1, for any temperature T, 
the dynamics of the system is essentially Markovian [59, 64], 
as the characteristic time of the system is much larger than 
the corresponding relaxation time of the bath. In contrast, 
the regime wo/wc > 1, where the characteristic time of the 
system is smaller or comparable to the bath relaxation time, 
is more akin to physically produce non-Markovian behavior. 
Moreover, we note that, as the evolution time tends to infin¬ 
ity, all non-Markovian signatures of the Gaussian channel are 
lost [59], irrespective of the system and bath parameters. Eor 
finite evolution times, the interplay between the diffusion and 
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Figure 6. (Color online) Non-Markovianity indicators in the quantum Brownian motion, plotted as a function of a, for different (scaled) 
reservoir temperatures: (a) T = 0, (b) F = 0.5, (c) T = I, and (d) T = 4. The non-Markovianity witness Nq is plotted for the squeezed 
thermal state with F 2 = 1 (red online dashed lines) and the mixed thermal state with ki = \ (blue online solid lines). All the states have fixed 
average number of excitations, h = 2.5. The Ohmic cut-off frequency and the characteristic frequency have been set at a>c = 1 and ojg = 6, 
respectively. The measure based on divisibility criteria, Afo, is plotted as well for comparison (black online dotted lines). All the quantities 
plotted are dimensionless. 


damping coefficients in determining the non-Markovian dy¬ 
namics, in the regime a>ofa>c > 1 , is dependent on the temper¬ 
ature T of the reservoir. In the high-7“ limit, A(f) » 7 (f), and 
the non-Markovianity of the evolution is solely determined by 
the diffusion coefficient; specifically by the condition A(f) < 0 
[56, 59, 64]. However, in the low-7“ regime, both A(f) and y(f) 
are comparable and the dynamics is intrinsically richer. In our 
investigation of non-Markovianity in QBM, we analyze both 
the low and high temperature regime, with the ratio ojq/oJc 
greater than unity. Furthermore, motivated by the analysis 
in the previous section, we only consider the non-Markovian 
witness Nq evaluated for the two classes of initial Gaussian 
states represented by Eqs. (27) and (28), with kj - k 2 - 1. 

In Figs. 5 and 6 , we have plotted Nq for values of the 

scaled temperature T (= ^^) ranging from 0 (vanishing) 
to 4 (high), with the ratio cjqIuc equal to 4 and 6 , respec¬ 
tively. We numerically calculated the non-Markovian witness 
for the two classes of initial Gaussian states, namely the en¬ 
tangled squeezed thermal state and the separable mixed ther¬ 
mal state. For comparison, we also plot in Figs. 5 and 6 the 
non-Markovian measure constructed from the violation of di¬ 
visibility, Ndi given by Eq. (15) [59]. The non-Markovian be¬ 
havior of the channel for QBM is evident in these figures, for 
both indicators No and Nq, obtained via non-divisibility and 
using GIF, respectively. We note that, apart from the case of 


small ratio wq/Wc at T = 0 case, the entangled two-mode STS 
with k 2 - I appears to give the maximum value of the non- 
Markovian witness among the considered initial probes. How¬ 
ever, the separable two-mode MTS with ki = 1 consistently 
provides a finite value of Nq in a broad parameter range, thus 
demonstrating once more that entanglement between the two 
modes in the input probe states is not necessary in principle to 
detect non-Markovianity of Gaussian channels, specifically in 
the case of the QBM. Recall that this was explicitly shown for 
the damping master equation in Sec. Ill D as well. 

Quite surprisingly, we observe that at T = 0, the witness 
Nq constructed from the entangled pure STS, in the regime 
- 4, does not show any perceptible value as compared 
to the one constructed from the separable MTS, or to the mea¬ 
sure Nd- Therefore, in this case, mixed separable states are 
necessary in some cases to detect non-Markovianity of zero- 
temperature QBM whereas entangled probes would fail to de¬ 
tect it, within our approach. A possible explanation for such 
a phenomenon may be gained by reaching a clearer under¬ 
standing of the interplay between the diffusion and damping 
coefficients in this regime. At very low T, up to first order cor¬ 
rections, the diffusion coefficient A(f) is independent of T and 
does not uniquely dominate non-Markovian characteristics in 
regimes close to resonance, wq ~ [56]. Hence, at T = 

0, the non-Markovian behavior is primarily dependent on the 
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value of 7 (f) with respect to A(f). This regime of imperceptible 
A(f) dependence, is not captured by the pure STS {k 2 - 1) in 
Fig. 5. However, as the ratio (OqIcjJc becomes higher and gets 
away from resonance, A(f) becomes negative and dominates 
the non-Markovian behavior, which is now duly captured by 
the pure STS in Fig. 6. 

In Fig. 7 we depict the interplay between damping and dif¬ 
fusion at r = 0, giving a pictorial exposition of the argu¬ 
ments presented above. We plot specifically the quantities 
(A(f) -H y(t))l2 and (A(f) - 7(f))/2, for 0 < toolcjc < 8. 
The figures show that the quantity (A(f) + 7(f))/2 is always 
positive at the regime closer to resonance, as compared to 
(A(f) - 7(f))/2 which features negative peaks. In this regime, 
the non-Markovianity is not dependent on A(f) [56]. However, 
in the regime oJoltOc » 0, both quantities (A(f) -i- y(t))l2 and 
(A(f) - 7(f))/2 can be negative, with A(f) having a stronger 
impact on the non-Markovianity. 


IV. CONCLUSION 

In this work we have presented a novel approach to 
witness—via Nq{A), Eq. (17)— and measure—via Nq{A), 
Eq. (18)— the non-Markovianity of a continuous variable 
Gaussian channel A, using bipartite Gaussian states as probes. 
Our framework is based on the breakdown of monotonicity of 
an operational figure of merit defined in the context of quan¬ 
tum metrology, the so-called Gaussian interferometric power 
[61, 62]. This work and Ref. [48], taken together, demon¬ 
strate that a general framework based on the interferometric 
power can be very fruitful to achieve an experimentally feasi¬ 
ble characterization and, in some relevant cases, an analytical 
quantification of non-Markovianity in dynamical maps span¬ 
ning from qubits to continuous variable systems. 

In the present paper we have applied our methods to study 
in detail two examples of open Gaussian dynamics, the damp¬ 
ing master equation and the quantum Brownian motion. 

Eor the damping master equation, we have calculated the 
measure Nq exactly, showing that pure entangled two-mode 
squeezed states are the optimal probes for the detection of 
non-Markovianity. Since our indicator is based on the quan¬ 
tum Eisher information, we may borrow a metrological termi¬ 
nology and argue that non-Markovianity can be detected at the 
Heisenberg limit, as reflected by the fact that Nq = ex¬ 

hibits a quadratic scaling with the mean energy of the probes. 
However, even separable states such as mixed thermal states 
are found to be useful to detect non-Markovianity, albeit re¬ 
sulting in suboptimal witnesses which scale at most 

linearly with the mean energy of the probes. 

For the quantum Brownian motion, we have analyzed dif¬ 
ferent regimes in terms of the reservoir temperature and 
the bare frequency of the system. We have shown that 
our approach successfully witnesses non-Markovianity for all 
regimes in which this behavior is expected to manifest, and 
our results are consistent with a recent characterization of 
non-Markovianity in terms of the non-divisibility criterion 
for Gaussian channels [59]. It is worth mentioning that we 
found some particular regimes (at zero temperature) where 
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Figure 7. (Color online) We plot the temporal evolution of the quanti¬ 
ties (a) (A(r) - y(t))l2 and (b) (A(f) -l- y{l))l2, under the master equa¬ 
tion for the quantum Brownian motion, in the T = 0 regime, as a 
function of time t and of the ratio (jJqI Notice how the quantity in 
(a) has several negative oscillations, while the quantity in (b) is only 
negative in a small region at short time and high frequency ratio. All 
the quantities plotted are dimensionless. 


non-Markovianity is better detected using initial states which 
are mixed and with no entanglement, as compared to pure 
two-mode squeezed probes. 

Overall, our approach demonstrates that cheap (in terms 
of engineering demands) two-mode resource states can be 
adopted to construct robust non-Markovianity witnesses for 
Gaussian channels acting on one mode of the system. 

A complete and practical characterization of non- 
Markovianity in continuous variable open systems, specifi¬ 
cally tailored to Gaussian states and Gaussian channels, may 
be of great importance for technological applications [43], es¬ 
pecially if one is interested in using memory effects and back- 
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flow of information from the environment, that can arise from 
non-Markovian evolutions, to protect and enhance quantum 
information tasks such as communication, cryptography and 
metrology. The role of non-Markovianity is also very promis¬ 
ing in quantum thermodynamics, since non-Markovian chan¬ 
nels can, in principle, allow work extraction accompanied by 
an increase of the system entropy [26]. We believe that our ap¬ 
proach can be further developed and applied to assess environ¬ 
mental enhancements in quantum technologies realized with 
continuous variable systems, as well as to reach a more funda¬ 
mental understanding of the physical mechanisms underpin¬ 
ning non-Markovianity in different realizations of Gaussian 
quantum channels. 
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